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ПРИВЕДЕНИЕ ВИХРЕВОГО МАГНИТНОГО ПОЛЯ 
К ПОТЕНЦИАЛЬНОМУ ПОЛЮ ИСТОЧНИКОВ 


В настоящей статье излагается методика расчета и моделирования стационар- 
ных и квазистационарных магнитных полей, основанная на преобразовании вихре- 
вого магнитного поля в потенциальное поле источников. В результате электриче- 
ские токи заменяются магнитными зарядами, а уравнения для векторов поля заме- 
няются уравнением для скалярной функции. 
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Введение. Метод приведения вихревого магнитного поля к потенциаль- 
ному полю источников основывается на тождественности магнитных полей 
системы токов и системы магнитных масс. 

Постановка задачи. Суть этого метода заключается в представлении ис- 


комого поля вектора Н в виде суммы расчетного поля Н, и вспомога- 
тельного поля Н.: 

Н=Н,+Но. (1) 
В результате введения вспомогательного поля становится возможным 


представить систему уравнений, описывающую стационарное магнитное 
поле, в виде двух систем 


тоНо= 8, (2) 
аи Но=-рт, 

ОЕ Н р = 0, (3) 
аи Нр=рт, 


одна из которых (система 3) является полным аналогом электростатическо- 
го поля, для которого можно использовать понятие скалярного магнитного 


потенциала. В дальнейшем задача определения поля вектора Н, заключа- 
ется в решении уравнения Пуассона: 

Фи ета4И „= -ри. (4) 

На выбор поля вектора Н, налагается единственное ограничение в 

виде условия хо Н, = 6. Многозначность этого уравнения может быть лик- 

видирована при задании закона распределения величины 4йжН». Одна- 

ко, если каким-либо образом выбрано удовлетворяющее уравнение 

коН,=$ поле Н,‚то значение величины Айя Н, в любой точке про- 


странства определено однозначно. Так, например, можно представить 
одно из решений, удовлетворяющее указанному условию в виде: 


фы __ 
Но= 9(6) а, (5) 
0 
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где (46- координатный вектор; @,- начальная линия отсчета от задан- 

ной точки, линии или поверхности; $(<6)- вектор плотности тока, 
являющийся в общем случае функцией текущей координаты @. 

При определении Н. по (5) возможен различный выбор направле- 

ния пути интегрирования, так как @, и с‹/6 выбираются произвольно. 


Именно некоторая свобода в выборе Н, дает возможность найти опти- 
мальный вариант расчетной модели, позволяющий упростить решение по- 
ставленной задачи, так как выбором поля Н, определяется объем вычис- 


лительной работы при нахождении поля Н,. Наиболее целесообразным 
направлением пути интегрирования является такой путь, который обеспе- 
чивает в конечном итоге отсутствие объемных магнитных зарядов р» и 


приводит лишь к возникновению поверхностных магнитных зарядов 0». 


Постановка задачи в таком виде сформулирована как первичная задача 
минимизации источников, регламентирующая выбор координатных по- 
верхностей так, чтобы вектор Н, был на них по возможности постоянным 
и имел составляющие, касательные к этим поверхностям. Так, например, 
при анализе существующих конструкций обмоток МГД- генератора целесо- 
образно в качестве исходной системы координат выбрать либо прямо- 
угольную, либо цилиндрическую систему координат. Однако варианты вы- 
полнения прямолинейных частей седлообразных обмоток указывают на 
широкое применение конструкций, имеющих цилиндрическую форму. 
Поэтому для дальнейшего анализа в качестве основной выбрана седлооб- 
разная обмотка, у которой ступенчатое распределение секций обусловлено 
как технологическими соображениями изготовления и сборки узлов магнит- 
ной системы МГД - генератора, так и для уменьшения максимального зна- 
чения магнитной индукции в лобовой части обмотки. При выбранной си- 
стеме координат расчет поля Н, не представляет трудности. Наибольший 


объем работ заключается в вычислении скалярного магнитного потенциала 


(/ 


777 , решение для которого в случае однородной среды может 


быть представлено в виде: 








1 р „ау б „45 р» созу 45 
[94 = т + т + т 
"4 з т г =? 1, (6) 


5; и 

где р„- момент двойного слоя магнитных зарядов; \ - угол между р» 

И р,а » направлен от точки приложения Р»„ к точке, в которой 
определяется С. 

При решении задачи минимизации областей с объемными магнит- 

ными зарядами р „, , на поверхностях, ограничивающих токовые зоны, воз- 

никает скачок напряженности магнитного поля. Физически это соответству- 
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ет случаю расположения на поверхностях простого слоя магнитных заря- 
ДОВ: 

0 т г Ри а 7: Но) (7) 
и поверхностного тока 

= ВоКН,-Н„,) , (8) 


где Н„,Н„,,Н„,,Н о, - нормальные и касательные составляющие 


вспомогательного поля вектора т на границах тОковой ЗОНЫ. 


Так как в реальной задаче на поверхности проводника источники 
отсутствуют, они Должны быть компенсированы магнитным слоем с плот- 
НОСТЬЮ 


‘„=-б„ (9) 

и поверхностным током 
К =-К. (10) 
Согласно правилу правого винта вектор плотности данного тока 
должен составлять угол 90° с направлением, касательным к вектору Н,. 
В целом ряде практических случаев для успешного решения задачи мини- 
мизации источников потенциального поля оказывается целесообразным 
введение двойных слоев магнитных зарядов, эквивалентирующих токовые 


слои. Такой искусственный прием обеспечивает ограничение эпюры Н„ и 
в значительной степени облегчает задачу исследования. Величина момента 
такого двойного слоя р„ определяется из условия, что скачок касательной 


составляющей напряженности при переходе через ДВОЙНОЙ СЛОЙ определя- 
ется скоростью изменения момента двойного слоя в точке перехода: 


Г др 
НЕ, = 
т ©; Г] 9. ' (11) 





где Р„- пр, Не, И На, - значение касательных составляющих вектора 


напряженности Н, соответственно с внутренней и внешней стороны 
двойного слоя. 
Отсюда 


Ри = НН». Ю Н»,) ав, . (12) 


и 
Решение уравнения (6) осуществляется с учетом найденного рас- 
пределения в пространстве простых магнитных слоев с поверхностной 


плотностью б„ и двойных магнитных слоев с моментом рР„, либо поверх- 
ностного тока /у. 

Каждый частный случай практического приложения данного метода 
отличается своими особенностями, обусловленными прежде всего реаль- 
ной геометрией системы и желанием упростить решение поставленной за- 
дачи. 
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Специфической особенностью рассматриваемой задачи является по- 
строение расчетной модели источников потенциального поля, которая 
позволяет одновременно достаточно полно учитывать реальную геометрию 
исследуемой седлообразной обмотки с учетом переходной области и в мак- 
симальной степени минимизировать количество этих источников. Необхо- 
димость последнего продиктована прежде всего требованием создания эф- 
фективной оптимизационной программы, в которой вопрос быстродействия 
расчета магнитного поля является одним из основных. Поэтому выбор оп- 
тимального варианта расчетной модели, из совокупности построенных, 
оказывается принципиальным при его численной реализации. 

Выводы. Алгоритм решения задачи расчёта магнитных полей при помощи 
фиктивных магнитных зарядов представляется следующим образом: 


» выбирается поле Н., которое удовлетворяет условию 
го Н,=0; 

® рассчитывается значение магнитных зарядов р» = - амМинН.; 

® решается уравнение @/мидгааИм = — б" для скалярного магнит- 
ного потенциала; 

® определяется поле Н› как Н› =- дгадИм и затем  иско- 
мое поле Н =Н, + Н,. 

В результате электрические токи плотностью к заменяются маг- 


нитными зарядами с плотностями р», о, ; уравнения для векторов поля за- 
меняются уравнением для скалярной функции Им, а вычисление искомого 


поля Н связано со сложением полей Н, и Н.. 
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0.3.РАВ$НТМ, М.В.УТМОКУВОУ 


ТВАМ$РОВМАТТОМ ОЕ А УОВТТСАЕ МАСМЕТТС ЕТЕГО 
ТО А РОТЕМТТАЕ ЕТЕГО ОЕ $ЗОЧВСЕ$ 


п ргезепЕ а|ае {Ле дезюп ргоседиге ап поде!пд оЁ {Пе $%аЧопагу апа а|- 
ппо5Ё $(афНопагу тадпейс Не, Базе оп {гапогтаНоп оР а уегса! тадпейс 
Не © а роепНа! оЁ зоигсеб 15 %ае4. Тп гезик, @есис сиггет аге герасеа 
МИ тадпеЧс спагдез, апа Ме едиаНопз Гог месбюг$ оЁ а Не аге герасеа 
М/И {Пе едцайоп Гог $са!аг Гипсйоп. 
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